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Despite the interest in Dark Matter (DM) searches is currently more focused on underground 
experiments, a signature of DM annihilation/decay in gamma-rays from the space would constitute 
a smoking gun for its identification. Here, we present the results of the survey of Segue 1 by the 
MAGIC-I telescope performed in 2008 and 2009. This source is considered by many as the most DM 
dominated Milky Way satellite galaxy known so far. The nearly 43 hours of data taken constitute 
the deepest observation ever made on a single dwarf galaxy by Cherenkov telescopes. No significant 
gamma-ray emission was found above an energy threshold of 100 GeV. Integral upper limits on 
the gamma-ray flux were calculated assuming various power-law spectra for the possible emission 
spectrum and for different energy thresholds. We also discuss a novel analysis that fully takes into 
account the spectral features of the gamma-ray spectrum of specific DM models in a SuperSymmetric 
scenario. 



INTRODUCTION 



In the ACDM cosmological scenario about 80% of 
the matter of the Universe is believed to be composed 
of non-baryonic matter, called Dark Matter (DM). 
The most popular DM candidates are the WIMPs 
(weakly interacting massive particles) supposed to 
be cold, electrically neutral, stable, and massive 
Among the huge plethora of WIMP candidates, the 
best motivated ones are related to the SuperSymmet- 
rical (SUSY) and Extra Dimensional extensions of the 
Standard Model of particle physics 0j- 
In the Minimal SuperSymmetric extension of the 
Standard Model (MSSM), the neutralino x represents 
an excellent cold DM candidate with a relic density 
compatible with the WMAP bounds. Since the neu- 
tralino is a Majorana particle, pairs of x can annihilate 
into Standard Model particles, e.g., quarks, leptons, 
and W bosons. The subsequent hadronization of those 
particles results in a continuum emission of gamma- 
rays characterized by a cut-off at the neutralino mass 
and by possible spectral features like bumps or a hard- 
ening of the spectral slope. 

The expected gamma-ray flux from DM-annihilating 



astrophysical objects, as function of the energy thresh- 
old Eq and the integration region Af2, within which 
the signal is integrated, can be factorized in two terms: 



$(> E ,An) = $^(> Eo)J(An) 



(1) 



The so-called particle physics factor $ pp depends on 
the features of the DM particle, and can be written 
as: 



<£ PP (> E ) = 



1 < <J„ 



4tt 
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(2) 

where < a ann v > is the velocity averaged annihilation 
cross-section, and B l is the particular branching ratio 
for the z-th annihilation channel. 
The term J (AO) (the so-called astrophysical factor) is 
given by the line-of-sight integral over the DM density 
squared within a solid angle Af2, and depends on the 
density profile of the DM halo of the source: 



J(Afi) = / / 
J An Jios 



p 2 (r(s,n))dsdn. 



(3) 



Since the gamma-ray flux of DM annihilation is pro- 
portional to the square of the DM density, only sources 



2 



201 1 Fermi Symposium, Roma., May. 9-12 



with high expected DM densities are good targets 
for DM indirect searches. Among these, the dwarf 
spheroidal satellite galaxies (dSphs) of the Milky Way 
(MW) are interesting objects thanks to their relative 
proximity to the Earth, to their high mass-to-light ra- 
tio (with values within tens and thousands of Mq/Lq) 
and to the expected absence of conventional gamma- 
ray sources within the system [3j, |4j. So far, around 
two dozen dSphs have been identified. Segue 1, discov- 
ered in 2006 by the SDSS 0, is located at 28 kpc from 
the Galactic Center, at (RA,DEC) = (10.12'\ 16.08°). 
Kinematics studies applied to 66 member stars al- 
lowed to estimate its mass-to-light ratio to be in the 
range 1320-3400 M Q /L @, highlighting Segue 1 as 
the most DM dominated dSph known so far. 
The MAGIC-I telescope is a 17 m dish Imaging Atmo- 
spheric Cherenkov Telescope (IACT), located at the 
Roque de los Muchachos Observatory, in the Canary 
Island of La Palma (2200 m a.s.l.). Thanks to its 
low energy threshold (~60 GeV at Zenith), high flux 
sensitivity (1.8% of the Crab Nebula flux in 50 hour 
of observations above ~250 GeV), and good angular 
and energy resolution (0.1° and 30% respectively, at 
100 GeV) Q, MAGIC-I is a suited instrument for the 
indirect search for DM candidates with energy of the 
neutralino mass or the Kaluza-Klcin state. 



II. MAGIC-I OBSERVATION AND ANALYSIS 
RESULTS 

A search for a possible DM gamma-ray signal com- 
ing from Segue 1 was performed by the MAGIC-I tele- 
scope between November 2008 and March 2009, for a 
total of 29.4 hours of observation time (after data se- 
lection). The data analysis was performed using the 
standard MAGIC-I analysis and reconstruction soft- 
ware [H . The number of gamma-ray candidates events 
from the direction of the source was estimated using 
the distribution of \a\ angles, which are related to 
the orientation of the showers. The overall analysis 
cuts were optimized and cross-checked for point-like 
sources with the aid of contemporaneous Crab Neb- 
ula data. In Figure [TJ the |a|-plot above 100 GeV is 
shown. The number of excess events was computed 
in a fixed fiducial signal region with \a\ < 14° and 
resulted to be N exc (> 100 GeV) = -279 ± 329, corre- 
sponding to a significance of — 0.85ct, computed using 
eq.(17) of Li&Ma 9]. Since results were consistent 
with no signal over the background, we derived Upper 
Limits (ULs) on the flux, calculated using the Rolkc 
method [HI at 95% confidence level, and assuming 
a 30% systematic uncertainty. Figure [2] shows the 
integral ULs achieved by the MAGIC-I observation 
of Segue 1 considering different energy thresholds Eq 
and different power-law spectra with spectral index 
r = -1, -1.5, -1.8, -2, -2.2, -2.4. It is worth noting 
that, using the Rolke method, the ULs on the num- 
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FIG. 1: a-plot from 29.4 hours of Segue 1 observation 
above 100 GeV. Red points represent the signal (ON distri- 
bution), black points the background (OFF distribution), 
and green points their difference. The vertical dashed line 
at a = 14° is the fiducial region below which the excess 
event number is estimated. Image taken from (TFJ]. 
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FIG. 2: Integral flux ULs from Segue 1. The arrows indi- 
cate the integral flux upper limits for different power-law 
spectra and energy thresholds. The dashed lines indicate 
the corresponding integral ULs if zero significance OLi,Ma 
is assumed. Image taken from [Iol |. 



ber of the excess events, and consequently the integral 
flux ULs, are affected by statistical fluctuations quan- 
tified by the significance of the observation a Li, Ma- 
This is an intrinsic feature of the statistical method 
exploited in the analysis and it should be taken into 
account when comparing ULs from different analyses. 
To show this effect, in Figured] we plot also the ULs 
(dashed lines) computed assuming a value for o~ Li. Ma 
equal to zero (with number of ON events equal to 
the number of OFF events in the signal region of |a|- 
plot) for different values of spectral index and energy 
threshold. 
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III. CONSTRAINTS ON DARK MATTER 
MODELS 

Assuming a particular form for Segue 1 DM halo, 
and a given particle model for the DM candidates, 
we can translate the integral ULs derived from the 
Segue 1 observation into constraints on the DM anni- 
hilation rate. 

Motivated by results from cosmological simula- 
tion, the DM halo around Segue 1 was mod- 
eled by using the Einasto radial profile [l2j 
with a s =1.1 xlO 8 Mq kpc -3 , r s =0.15 kpc, and 
n=3.3. With those parameters, the total astro- 
physical factor of Segue 1 results to be J(Ail) = 
1.78xl0 19 GeV 2 cm~ 5 sr. Since the analysis was per- 
formed assuming point-like source cuts (corresponding 
to an angular integration of 0.14° above 100 GeV), 
we estimated the effective astrophysical factor within 
the analysis cuts to be used in the following analy- 
sis, being its value J(AO) = 1.14xl0 19 GeV 2 cur 5 sr 
(corresponding to the 64% of the total astrophysical 
factor). 

Concerning the particle physics, we restricted our- 
selves to the case of a SUSY model in which the pres- 
ence of a discrete symmetry (R-parity) guarantees 
that the Lightest SuperSymmetric Particle (LSP) is 
stable over cosmological timescales and, therefore, a 
good DM candidate. We considered a 5-dimensional 
subspace of the MSSM called mSUGRA [l]|, for which 
the basic parameters are the universal masses of the 
gauginos (m^) and scalars (too), the trilinear cou- 
pling (Ao), the ratio of the vacuum expectation val- 
ues of the two Higgs fields (tan/3) and the sign of the 
Higgsino mass term (sign(/i)). In order to study the 
phenomenology of mSUGRA we performed a grid scan 
over the parameter space, for a total of 5xl0 6 points 
(for the details see |10|). 

The full circles of Figure [3] represent all the models 
of the scan i) where the lightest SUSY particle is a 
neutralino, ii) that survive the Standard Model con- 
straints and Hi) with a relic density compatible with 
the value derived by WMAP data within three times 
its experimental error aw map [l4{. For each DM 
model of the scan, we computed the integral flux UL 
(above an energy threshold E ), using the Segue 1 
data and the specific gamma-ray spectrum derived 
from the individual DM model. Since the spectra for 
each DM model have different shapes and cut-offs, the 
value of the optimal energy threshold Eo was com- 
puted individually for each DM mass. We then con- 
verted the flux ULs into ULs on the velocity averaged 
cross-section to have a direct comparison of exper- 
imental data with the theoretical predictions. The 
results are plotted in Figure [3] as function of the neu- 
tralino mass: each DM model of the scan (full circle) 
is compared to its own UL (square). For each point 
we defined an enhancement factor (ENF) as the ratio 
between the UL on the velocity averaged cross-section 



Upper limits 
mSUGRA (<WMAP) 
mSUGRA (WMAP) 
ENF<10 4 , 
10"<ENF<1Cf 
10 <ENF<10 
ENF>10 




MAGIC Coll. (g011) 



800 1000 1200 

m, [GeV] 



FIG. 3: Velocity averaged annihilation cross-section ULs 
from Segue 1 MAGIC-I data computed for individual 
points in the scan. Grey crosses indicate the velocity aver- 
aged annihilation cross-section value for those points in the 
scan that pass the SM constraints and with a relic density 
lower than WMAP bound. The full circles only consider 
models within Saw map from WMAP bounds. For each 
of these full circles the UL on the cross-section can be 
computed from the Segue 1 data (after energy threshold 
optimization) and it is indicated here by a square. Cir- 
cles and squared are colored in term of the enhancement 
factor. Image taken from [To| ] . 



and the value predicted by mSUGRA. This quantity 
quantifies how far away we are from excluding some 
portions of the mSUGRA parameters space. From 
Figure [3] it can be seen that ENFs for model compat- 
ible with the WMAP bounds are typically above 10 3 , 
while typical values are of the order of 10 4-5 . 



IV. IMPACT ON PAMELA PREFERRED 
REGION 

We tested our ULs on some of the models pro- 
posed in the literature that can explain the PAMELA 
data [l5| for the energy spectrum of the positron frac- 
tion e + /(e + + e~) as due to DM annihilation into 
leptons. The regions in the (m x , < a ann v >) plane 
that provide a good fit to the PAMELA measure- 
ments [TH, [l?} for three different channels of DM 
annihilation are shown in Figure HI The annihila- 

t + t~ have been 



XX 



tion channels XX A <+ M _ j 
taken from SuperSymmetry, while for \X ~* 
2e + e~ the existence of a new dark force, mediated by 
the carrier <fi that decays into leptons [lij, has been 
assumed. In Figure [4] we plot the ULs obtained from 
the Segue 1 data, using again the specific DM anni- 
hilation spectra. We can see that, in this case, the 
ENFs needed to meet the PAMELA-favoured region 
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FIG. 4: Exclusion lines for a neutralino DM annihilating 
exclusively into j-L + )J,~ (green lines) or t + t~ (blue line), 
and for a DM candidate interacting with a light interme- 
diate state (f> decaying into a pair of electrons (pink line). 
The same annihilation channels (with the same color cod- 
ing) are considered to draw the regions in the plane that 
provide a good fit to the PAMELA measurement of the en- 
ergy spectrum of the positron fraction. Image taken from 



are much smaller than those found for mSUGRA, and 
in the case of annihilation into t + t~ our ULs are prob- 
ing the relevant regions. However, it is worth mention- 
ing that, since the uncertainty in the Segue 1 astro- 
physical factor is quite large [12], an improvement in 
accuracy of astrophysical factor value could be able 
to put more stringent constrains and to confirm the 
exclusion of the PAMELA region for DM particle an- 
nihilating in t + t~. 

V. CONCLUSIONS 

A search for a possible DM gamma-ray signal com- 
ing from Segue 1 was performed by the MAGIC-I tele- 
scope. No hints of signal were found above the back- 
ground for energies larger than 100 GeV. Integral ULs 



on the gamma-ray emission were computed assum- 
ing different power-law energy spectra. Within the 
mSUGRA scenario, a large scan of neutralino mod- 
els was performed over the parameter space. Sub- 
sequently for each simulated DM model, the ULs 
on the velocity averaged annihilation cross-section 
(m x ,< a ann v >) were derived separately for each 
point in the scan in order to account for the depen- 
dence on the specific spectra. Results indicate that a 
general exclusion plot cannot be drawn to constrain 
the parameter space, so we provide the results in terms 
of enhancement factors. A minimum boost on the flux 
is found of the order of 10 3 (for models compatible 
with WMAP) while the typical values are at 10 4 " 5 . 
MAGIC-I data of Segue 1 can be useful to put con- 
straints on those DM models that are provided in the 
literature to explain the PAMELA data. Our ULs are 
probing the PAMELA region for the DM models an- 
nihilating into t + t~ but the robustness of this result 
could be improved, decreasing the uncertainty in the 
astrophysical factor. 

Although the MAGIC-I observation did not result in a 
detection, and the ULs require still high flux enhance- 
ment factors to actually match the experiment sensi- 
tivity, an analysis like the one presented here is able to 
point out details and features that can be important 
for future deep exposures of this or similar objects, 
with next-generation Cherenkov experiments. 
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